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Abstract 

We examine extended bound states in the dilute limit of the extended Hubbard 
model with next-nearest-neighbour (nnn) hopping and nnn interaction on the 
two dimensional square lattice. By the exact solution of a two-body problem we 
show that d-wave pairing is strongly favored by the nnn hopping. We determine 
the binding energies, mobilities and dispersion curves across Brillouin zone (B.z.) 
for bound states of various symmetries. The numerical analysis shows that the 
deep bound states can exist in the whole B.z. and we have also found s* — d 
mixing of the bound states. 



I. INTRODUCTION. 

The question of the symmetry of the superconducting order parameter arose very soon af- 
ter the discovery of high temperature superconductivity (HTS). Different symmetries imply 
superconductivity of considerably different properties and large efforts have been devoted to 
determining the type of pairing in different materials. Experimental findings have yielded am- 
biguous results, yet they strongly suggest explanations extending beyond the simplest s-wave 
picture, possibly (i-wave or mixed state in cuprate HTS. Recently, the angle-resolved photoe- 
mission spectroscopy and tri-cristal Josephson effect [2] revealed a large (i-wave component in 
the order parameter of cuprate superconductors. Apart from different pairing mechanisms fa- 
voring this kind of pairing, this situation may be due to kinetic mechanisms, as suggested by 
Blaer et al [1]. With inclusion of the nnn hopping, a bound pair formed on adjacent lattice 
sites can propagate over long distances before decaying. The effective mass of such a pair can 
become even smaller than the mass of its constituent fermions. The issue of nnn hopping has 
also important implications for the problem of hole motion in an antiferromagnetic background 
[3] and for the explanation of the HTS properties [4]- [6]. 

Exactly solvable model of two fermions on a periodic lattice would give a valuable insight into 
the problems of both: the role of nnn hopping and the symmetry of the superconducting state. 
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It is also of importance, as a limiting case for the two relevant models of superconductivity: 
BCS and preformed pairs scenario. 

In this paper we report results of exact calculations of bound states (BS) of fcrmions on 
a 2D square lattice including nnn hopping and nnn interactions, using formulation of Ref. [1]. 
We have evaluated the pair binding energies, pair mobilities and dispersion curves across the 
Brillouin zone for the extended Hubbard model for different pairing symmetries. We discuss in 
detail the effect of nnn hopping t2 and nnn interactions W2, on the stabihty of the shallow BS. 



II. FORMALISM 
A. The model 

The details of the theoretical treatment are given in [1] so we shall give here only the main 
points. We start with the generalized Hubbard Hamiltonian on a periodic lattice with N sites: 

H = Ho + H^, (1) 
Ho= E t{f-f')cUf')ca{r-), (2) 

r,f',(7 

Hi = E[9sif- f')S^ir, r')S{r, r') + 9Tif- f')T.Tl{r, fWr, f% (3) 

r,f' 

where r, r ' denote the positions of electrons and S, T are singlet and triplet pairing operators: 

g(^ff ^/-) ^ c^ir)ci{f') + c^{f')ci{f) 

To{f, r') = ^T(r1q(fO-c,(fOc,(r-) ^ 

T,(r,r')^'J^^, (6) 
= (7) 

In the standard Hubbard model notation gs{Ri,Ri) = U. gs{Ri,Rj) = Wij and gx{Ri, Rj) = 
Wij. The inclusion of antiferromagnetic exchange J2ij Jij ' Si - Sj, Jij > amounts to change 
gs{Ri,Rj) = Wij — 3Jij/2, gT{Ri,Rj) = Wij + Jij/2 and triplet pairing can be suppressed. 
Prom now on we will consider the singlet pairing only. 

A convenient basis set in the Hilbert space of two fermions consists of states \P, f, o"i, (72 > 
with total Bloch momentum P and relative separation r. In relative coordinates in the reciprocal 
space the eigenvectors of the kinetic part of the Hamiltonian (1) are the states: 

|P,p,<7l,a2 >= clp^2)-P,aAp/2)+pJ^ >' 

with the eigenvalues: 
where 



ep- = —2ti{cos{pxS) + cos{pyS)) — 4^2 cos{pxS) * cos{pyS). (10) 

where 5 is a lattice constant. We also define 

1 



^{f\p) =< P, r, (71, a2\P,p, (71, (72 >= e'^•^ (11) 
^(r) =< P,r,(7i,(72|P,V,(^i,(72 >, (12) 



N 

1 



iV 

where ilj{r) is a wave function in the relative position space, a result of projection of any two- 
fermion state on a state with definite relative separation r and total Bloch momentum P. To 
solve the Schrodinger equation 

{Ho + H,)\P,'ilj>^E\P,'iP>, (13) 

where Hq is the kinetic part of the Hamiltonian (1), we introduce: 

|P,V;>=^F(p)|P,p>, (14) 
p 

which yields an integral equation for coefficients F{p): 

(15) 

{E - Ep.)F{p) = J2FiQ) < Pp\Hi\PQ>- 

q 

r q 

Using the above formula for F{j)) and definition: '^(r) — J2pF{p)i/j{f\p), one obtains the 
final equation for the wave function ilj{r): 

V'(r1 ^Y.G{E,p,r,r')g{r')i;{r'), (17) 

r' 

where: 

q Pq 

is the Green function. 

For the sites where interaction g is present, Eq. (17) can be cast into the form: 

[l-g{E,P)g]ij = 0, (19) 

where ip is a n x 1 matrix with ipj = ip[Rj), Q{E,P) is a n x n matrix with Qij{E,P) = 
G{E, P, Ri, Rj) and g is a diagonal nxn matrix with gu = g{Ri) where g{Ri) ^ only for sites 
r = 0, ±i?i, ±Rn-i. The energy eigenvalues are calculated from the determinant: 

det[l - G{E, P)g\ = 0. (20) 

To get the wave function for lattice sites outside described range we have to use Eq. (17). 



B. Final equations 



The matrix element in the reciprocal space is given as: 

N < Pp\H,\Pq>= go + 5i7p^-V + , (21) 

where: — 2[cos{pxS)+cos(py5)], 7?^ = 4:COs{px5) cos(py5), go — U is the on-site interaction, 
gi — Wi and g2 — W2 denote intersite interactions between nn and nnn, respectively. A 
separation of the matrix element (21) into parts depending exclusively on p and q is done 
by using five components ipnip), transforming according to irreducible representations of the 
symmetry group of the d — 2 lattice C^^ (given in the third column) : 

ipo = 1, (s-wave) Ai 

ipi{p) = cos{pxS) + cos{pyS), (ext. s*-wavc) Ai 

'4'2{p) = cos{pJ) - cos{pyS), (42_j^2-wave) Bi (22) 

ip3{p) = 2cos{pJ) * cos{py5), {sxy*-wecve) Ai 

ip^ip) — 2sin{pxS) * sin(pyS). (dxy-wave) B2 

With such a form of the basis functions the Green's function is given by: 

Above decomposition allows to write the following equation for F{p) being the Fourier transform 
of ■0(r)- 

F^^J29nGnm{E,P)F^, (24) 

for each irreducible representation, where = (1/-/V) J2pF{p)i'm{p)- Equation (20) for eigen- 
values E takes the form shown below - for s*-wave pairing (compare also [5,6]): 

l-Goof/ -GoiWi ~Go3W2 
det I -Goi^ 1 - GnWi -Gisl^a | = 0, (25) 
-GosU -GisM/i 1-G33W^2 



for ci3.2_2^2-wave: 



and for d^y'- 



l + |W^i|G22 = 0, (26) 



1 + 11^21^44 = 0. (27) 

The above decomposition holds for arbitrary t2 and P = or t2 = and P on the TM line 
as well as for any t2 with P on the TM line and W2 = 0. In general case one has to resort to 
Eq.(20). 

Eq. (25) reads: 



G^Goo - G^iGaa) - Gl^UW^ - Gl^W^W^ - Gl^UW^ = 0. (28) 



From Eq. (28) we can derive simpler cases [6]: 

Gn = -l/|iyi|, forW^2 = = C/; 

-(700(1 + GiilM^il) + G^ilW^il = 0, for = 0, C/ = oo. 

The relative pair mobility which is defined as: 



(29) 



(30) 



where m* is an efi^ective mass of a bound pair, and for a = ^2/^1 7^ 0, for two-dimensional 
square lattice 2m/ = h^/ih + 2t2W. Effective mass m* can be obtained by expanding the 
denominator in Eq.(23). We take: 

E - Ep^ = 2eo - Eb-ujp- ep/^-p - > (31) 

— * 

where E^ is the binding energy at P = 0, o;^ is the dispersion of the bound state and ep is given 
by Eq. (10). In the long-wavelength limit on the FM line: up = DP^ where D = h^S'^/2m,*. 
The coefficient D is obtained by expanding ep^2±p i^^o powers of P, inverting the series and 

substituting resulting Gnm{E,P) into appropriate Eqs (26 - 29). The final equations for pair 
mobilities in the simplest cases are: 

1 A 

(32) 



1 + 2a B, 



nn 



where n = 1 for s*-wave pairing with U = and W2 = 0; n = 2 for d^-i-y^ and n = 4 for dxy 
pairing with a — 0, respectively. We have used denotations: 

o _ J_ \- V'n(p)t^m(p)V^(p) .V 

A^^(2eo-i^;,-2e,)2' ^""^^ 

Bnm — BnmO- (35) 

Mobilities in other cases have more complicated forms. 



III. RESULTS 

In Figs (1) - (3) we have plotted the binding energies and mobilities vs |W^i/8ti| for extended- 
s pairing (denoted as s*) and for d3.2_j^2-wave pairing (in short d), for nnn hopping with param- 
eter a — a — —0.2 and a = —0.45. The last two values roughly correspond to the values of 
the parameter a, suggested by the band structure calculations in La2-xSrxCu04 {a ~ —0.25) 
and FBaaCMgOr-^ (a ~ -0.45) [4]. 

With only the nn hopping taken into account (i.e. a = 0) for U = infinitesimally small 
value of iH^il can create a bound state (of the s* type) - there is no critical value. For finite U or 



W2, including the case of the infinite on-site repulsion, \Wi\ has to reach critical value to create 
such a bound state (|VFic|/8ti = 0.25 for s*-wave for [/ = oo for a — 0). The d-wave pairing 
does not depend on U, and also a critical value of |Wic|/8ti = 0.915 (or |Vr2c|/8ti = 0.8268 in 
the case of d^y) has to be reached to form a pair of the d type. In this case (a = 0), the binding 
energies of d type are smaller than the s-wave binding energies in the whole range of displayed 
parameters. The mobility of pairs is smaller then 1 (i.e. the effective mass of a pair is larger 
then the mass of its constituent fermions) and decreases with increasing intersite attraction. 
Increasing \Wi\ makes the binding stronger, but at the same time, the pair movement more 
difficult. The mobility increases with decreasing \ Wi \ for a = 0, and is always smaller for rf-wave 
than s-wavc. Let us note, that the mobility of the s*-pairs for [/ = cxo becomes larger than the 
mobility of the s*-pairs for ?7 = near the critical value of Wi. It is so, despite the fact, that 
the binding energies for [/ = oo are smaller than in the U — case. It is also of interest, that 
both mobility and binding energies of the s*-pairs for U — oo = W2 resemble the respective 
quantities of the d-wscve pairing. Differences appear near the critical value, which is slightly 
smaller in the s* case, with s* mobility there about twice as large as in the d case, but they 
diminish very quickly with increasing \Wi\. 

The situation changes rapidly with increasing t2 (which is taken with the opposite sign than 
ti). The binding energy of ci-wave pairs grows and can become larger than that of s-wave pairs 
for large enough \ Wi\. The same happens for the mobilities. For the d-wave, the mobility grows 
faster than s-wave and for large enough \ Wi\ (although not the same as for pair binding energies) 
the mobility of (i-pairs exceeds that of s-pairs. Starting from \a\ ~ 0.3 the pair mobilities can 
become greater than 1 and even strongly bound pair can easily move without breaking its bond. 
For \Wi \ — > 00, Z ^ ±|a|/(l -|- 2a) with -|-(— ) for d (s*) -wave pairs, respectively [1]. 

The increase of pair mobilities with the parameter \a\ is much more abrupt for rf-wave 
pairing than for s-wave, so at the value a = —0.45 the maximum of rf-pairs mobility is around 
6 while that of s-pairs around 4, and the mobility for s-wave is smaller than for (i-wave in the 
whole range of displayed Wi values (Fig. 3). 

With growing |a| a new feature develops: a maximum appears in the mobility curves. The 
mobility generally increases with a and decreases with \ Wi \ but is not the largest for the critical 
value of Wi. For s*-wave pairing with U = the mobility for Wi = is always equal to the 
mass of two free fermions, i.e. 1 in our units. For small binding energies, the mobility increases 
with \Wi\ and diminishes only later, for still larger \Wi\, which creates a maximum. This is 
a consequence of the saddle points located at (0, tt) and (tt, 0) in the reciprocal space. The 
dispersion curves are anomalously flat in that region (as can be seen in Fig. 5), particles move 
slower and the probability of pairing increases. The singularity is only "felt" by the shallow 
states, which results in a maximum for small \Wi\. 

We note that for large enough \ Wi\ the mobility of s* bound pairs becomes negative, which 
is connected to the fact that the minimum of BS energy moves from k — towards k — tt. 

In Figs (4) and (5) we show the bound state dispersion curves across the Brillouin zone for 
W2 = 0. The BS energies of d- and s*- pairing are the same at the zone boundary. At the F 
point d-wave BS energy is larger or smaller than s-wave BS energy, depending on a and \Wi\. 
The d-wecve BS energies can be larger than s-wave ones in a part of the B. z. and smaller in 
another part of the B. z. and the possibifity of mixed s* — d state emerges. 

In Fig. (5) we plotted the dispersion curves along the lines F — Af — X — F for a = —0.45. 
We observe a mixing of s* and d wave states and well developed dispersion for deep bound 



states. Calculations for small binding energies show, however, that shallow bound states may 
exist only in parts of the B.z. We were unable to to find bound states in the domain of real 
energies near the M point for small \Wi\ (both in the case of s*- and d-wave pairing). There 
is a possibility that BS enter the scattering continuum and turn into resonance states. The 
problem of the resonance states is left for further study. 

Fig. (6) shows the ground state diagram for the shallow bound states. The nnn interaction 
W2 was assumed to be positive, i.e. we disregarded d^y pairing. The resulting curves show, 
that for nn hopping s*-wave pairing prevails even for the infinite nnn repulsion. The inclusion 
of t2 together with the on-site repulsion act in favour of rf-wave pairing. For a given U , a finite 
W2 is needed to get the ground state of rf-wavc symmetry. 

Finally, let us point out that our results for existence of two-body BS in the s-wave channel 
are directly related to the stability of s-wave superconductivity on 2D lattice in a dilute limit [7] . 

IV. CONCLUSIONS 

We have examined two-electron bound states in an extended Hubbard model with nnn 
hopping and nnn interaction. The analysis shows that t2 affects kinetic energy of pairs in a 
considerable way. This second hopping, with opposite sign to ti, enhances (i-wave pairing more 
than s*-wave. The nnn hopping increases the relative pair mobility, which can become greater 
than that of free fermions (and a possibility of long-lived higly mobile bound pairs arises), in 
agreement with [1]. The increase of (ia;2_j^2-wave pair mobility is much larger than the respective 
one of the s*-wave. The nnn interaction influences the d-wave pairing only of the d^-y type. We 
have also exemplified the possibility of existence of a s* — dy.2_y2 mixed state. We have also 
pointed out that for large t2 shallow bound states may exist only in certain regions of B.z., in 
contrast to the deep bound states which can occur in the whole B.z., for strong enough intersite 
attraction. 

V. ACKNOWLEDGMENTS 

This paper was supported by the State Committee for Scientific Research (KBN Poland): 
Project No 2 P03B 056 14. 



REFERENCES 



[1] A. S. Blaer, H. C. Ren and O. Tchcrnyshyov, Phys. Rev. B 55, 6035 (1997). 

[2] Z. X. Shen et al., Phys. Rev. Lett. 70, 1553 (1993); C. C. Tsuei et al., Phys. Rev. Lett. 

73, 593 (1994). 
[3] A. Shiller et al, Phys. Rev. B 51, 8337 (1995). 

[4] P. Benard, L. Chen and A. - M. S. Tremblay, Phys. Rev. B 47, 15217 (1993). 
[5] A. Teubel, E. KoUey and W. Kollcy, phys. stat. sol. (b) 157, 389 (1990). 
[6] R. Micnas, J. Ranniger and S. Robaszkiewicz, Rev. Mod. Phys. 62, 112 (1990). 
[7] M. Randeira, J.-M. Duan, L.-Y. Shieh, Phys. Rev. B 41 327 (1990). 



FIGURES 
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FIG. 1. Binding energies (curves increasing with |II^i|) and pair mobilities (curves decreasing with 
\Wi\) on d = 2 square lattice with nn hopping (a = 0) for s*-wave pairing: for U = 0, W2 = 0; for 
U = 0, W2 = 00; for U = 00, W2 = 0; for U = 00, W2 = 00; and for d-wave. 
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FIG. 2. Binding energies (curves increasing with \ Wi\) and pair mobilities (curves decreasing with 
\Wi\) on d = 2 square lattice with nnn hopping (a = —0.2) for s*-wave pairing: for U = 0, W2 = 0; 
for U = 0, W2 = 00; for U = 00, W2 = 0; for U = 00, W2 = 00; and for cZ-wave. 



a=-0.45 




|WJ/8t, 

FIG. 3. Binding energies (inset) and pair mobilities for o = —0.45 for s*-wave pairing for 
U = = W2, for [/ = 0, W2 = 00, for U = 00, W2 = 0, for U = 00 = W2, and for d-wave. 
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FIG. 4. Dispersion curves for d-wave and s*-wave bound pairs for a = —0.2. Heavy line denotes 
the lower boundary of scattering continuum: min^ ^p0- 9 = -P is the total momentum. 
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FIG. 5. Dispersion curves for d-wave, s*-wave and mixed s* — d type bound pairs for a = —0.45 
along the lines F - M - X - T. 




